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' We discuss a quantum-metrology protocol designed to estimate a physical parameter in a Bose- 

Einstein condensate of atoms, and we show that the measurement uncertainty can decrease faster 
than 1 /N. The 1 /N scaling is usually thought to be the best possible in any measurement scheme. 
From the perspective of quantum information theory, we outline the main idea that leads to a mea- 
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P3 ' surement uncertainty that scales better than 1 /N. We examine in detail some potential problems and 

a ■ 

^ ' challenges that arise in implementing such a measurement protocol using a Bose-Einstein conden- 



sate. We discuss how some of these issues can be dealt with by using lower-dimensional condensates 



, trapped in nonharmonic potentials. 
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I. INTRODUCTION 

In quantum metrology, the description "Heisenberg-limited scaling" refers to the best possible scaling 
of the measurement uncertainty with the resources put into a measurement. The phrase arises not from 
Heisenberg uncertainty relations, but from uncertainty relations of the Mandelstam-Tamm type yj, 

dy{A^K}''^ > ^ , (1.1) 

in units with h = I. The uncertainty 6y in a parameter y that, in part, determines the state of a quantum 
system is related to the standard deviation of the operator K that generates translations of the state along a 
path parameterized by y. A sequence of logical and mathematical steps is needed to provide a rigorous con- 
nection between the problem of measurement uncertainty in quantum metrology and uncertainty relations 
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of the Mandelstam-Tamm type. The pioneering work of Helstrom (2^, Holevo ||3|], Braunstein, Caves, and 
Milbum ||4],|5|], and others laid out and elucidated these steps. We summarize them below for the sake of 
completeness. 

The discussion in this paper is restricted to single-parameter estimation. The first step in estimating the 
value of a parameter is to identify an elementary physical system that is sensitive to changes in the parameter, 
just as one would choose a balance to measure weight or a thermometer to measure temperature. One or 
more of these elementary systems make up the measuring device or probe. The measurement uncertainty 
is a property of this measuring device. In quantum metrology this means that we expect the measurement 
uncertainty to depend on the initial state of the quantum probe, its evolution, and the measurement made on 
the probe to extract information about the parameter. The quantum Cramer-Rao bound quantifies the idea 
that the optimal measurement uncertainty is inversely proportional to the change in the state of the probe 
corresponding to small changes in the value of the parameter: 

Here i/i'DO denotes a distance element in the space of density operators of the probe, and 3(7, t) is the 
quantum Fisher information. The uncertainty in determining 7 is quantified by the units-corrected, root- 
mean-square deviation of one's estimate of the parameter, 7est, from the true value 7: 

-yfr. (1.3) 



^|5(7est>/57l 

In classical statistics, the Cramer-Rao bound on measurement uncertainty is given by 

1 



(Syf > -— , (1.4) 



where 



J(r)^|(^lnM^l7))y (1-5) 

called the Fisher information, is an average over the probability distribution p{^\y) for a random variable ^ 
and is a measure of the information that ^ can provide about 7. The classical Cramer-Rao bound (11.41) can 
generally be achieved only asymptotically in a large number of trials, i.e., independent measurements of 
^. The requirement of many trials to achieve the Cramer-Rao bound is important, but as a purely classical 
effect, it is not germane to our discussion of quantum limits, so we do not consider it further in the remainder 
of this paper. 

If the probe used to estimate the value of 7 were a classical system, then ^ would label the possible states 
of the probe at the end of the measurement, with p{^\y) being the probabihty of finding the probe in each 



3 



of these states. For a quantum probe in a state p{y, t) at the end of the measurement process, ^ labels the 
possible outcomes of a measurement performed on the probe, which is described by POVM elements E{£), 
with J E{^) = 1 and p{^\y,t) = tr[E{^)p{y,t)]. The classical Fisher information I{y,t), defined using 
pi^ly, t), clearly depends on the choice of POVM. The quantum Fisher information, which is independent 
of the choice of POVM, is therefore defined as 



3(7, f) = max I{y, t) . 



(1.6) 



The maximization over allpossible measurements in the above equation is a rather daunting prospect, but 



it can be shown that IB yL 



3(7, t) = tr[p(7, t)2\y, t)] = {2\y, t)) . 



(1.7) 



The symmetric logarithmic derivative, £(7, t), is the Hermitian operator defined implicitly by the equation 

(1.8) 



We now make two simplifying assumptions. First, we assume that translations in the parameter are 
generated by a unitary operator. This allows us to characterize the translations in terms of a Hermitian 
generator K{y, t) defined by 

dp{y,t) 



dy 



= -i[K{y,t), piyj)] . 



(1.9) 



Second, we assume that the state of the probe is pure, which implies p^{y, t) - p{y, t). Under these two 
assumptions, we can identify the symmetric logarithmic derivative as 



2(y,t) = l^ = -li[K{y,t), p(7,0] , 



(1.10) 



and the quantum Fisher information becomes 



3(7, t) - AtripK^ -pKpK) = 4{A^ K{y, t)) . 



(1.11) 



Thus, using Eq. (11.21 ). we obtain 



6y > 



1 



(1.12) 



2<A2/:(7, 0)1/2 ' 

which is a rigorous statement of the Mandelstam-Tamm uncertainty relation dl.ll ). Our two simplifying 
assumptions can be relaxed IS yL |4j, but we do not need the more general forms of the Cramer-Rao 
bound in this paper. 
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We can further simplify Eq. (11.121 ) by noting that the variance of a Hermitian operator is bounded from 
above by (A^K) < \\K\\^/4, where || • || is defined as the difference between the largest and smallest eigen- 
values of a Hermitian operator (this is a semi-norm for Hermitian operators). The quantum Cramer-Rao 
bound then becomes 

5y> ^- . (1.13) 

We mentioned the elementary quantum systems, sensitive to y, that are used to build the probe. The 
number N of such elementary units of the probe can be regarded as the most significant resource that goes 
into a measurement scheme. The differences between the tensor-product state space of a composite quantum 
system of N probe units and the Cartesian-product state space of an equivalent, classical composite system 
is the motivation for investigating whether a composite quantum probe offers advantages over classical ones 
in the relationship between 5y and A'^. 

From Eq. (11.131 ) we see that theoretically the A'^-dependence of the bound on 6y comes solely from the 
dependence of the generator K on A'^. In any particular quantum-metrology scheme, however, this bound 
might not be achievable, and additional dependence of 6y on N can come from the nature of the state of 
the probe as well. To see these dependences clearly and to understand what "Heisenberg-hmited scaling" 
means, we view quantum metrology from the perspective of quantum information theory using the language 
of quantum circuits in Sec. |lll We also explain how one can construct measurement protocols in which dy 
scales with N ina. manner not thought to be possible until recently. Sectionlllllexamines in some detail how 
such an enhanced metrology protocol might be implemented in a Bose-Einstein condensate (BEC) of 
atoms and considers the various problems and issues that might arise in a BEC realization of the proposed 
metrology scheme. 

II. QUANTUM METROLOGY FROM AN INFORMATION-THEORETIC PERSPECTIVE 

In this section, we follow Giovannetti, Lloyd, and Maccone |l6[] in using quantum circuits to describe 
and analyze metrology protocols. From this perspective, we first look at a couple of well-known measure- 
ment schemes that were considered in ^ — ^Ramsey interferometry (Sec. Ill Al ) and interferometry using a 
Schrodinger-cat state (superposition of macroscopically distinct states) (Sec. IIIBI ) — with the aim of gener- 
alizing these circuits to new protocols that were introduced in [7, 8]. In these initial discussions of Ramsey 
interferometry and cat-state interferometry, we assume that the elementary quantum systems that make up 
the probe are qubits. The quantization axis is taken to be along the z-direction of a Bloch-sphere repre- 
sentation, with the standard basis states along this direction denoted as |0) and |1). Despite the notation 
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FIG. 1 : Quantum circuit for Ramsey interferometry. 

the actual qubits need not be spin- 1/2 paiticles; they could very well be atoms in which only two energy 
levels are relevant or a variety of other suitable systems. In the subsequent general discussions of linear and 
nonlinear interferometry (Sees. HTCl and III Dl) . we allow the probe units to be any quantum system. It turns 
out, however, that optimal sensitivities are always attained by using only two levels of each unit, so in the 
end we can always regard the probe units as qubits. 

In all the quantum circuits depicted in this section, we use N = 3 probe units as an example. 



A. Ramsey interferometry 

A typical Ramsey interferometer, such as the one in |l9|], can be represented by the quantum circuit in 
Fig-El In this measurement protocol, each of the qubits that make up the probe evolves independently. All 
the qubits are initialized in the state |0), which might represent the ground state in Ramsey interferometry 
using atoms. The Hadamard gate H puts each of the qubits in an equal superposition of the two basis states, 
(|0) + U))/ V2. The parameter-dependent evolution of the quantum probe is generated by the Hamiltonian 

N 

^Ramsey ^ y^^^z; j/^ = jJz , (2-1) 
7=1 

where cr^. y denotes the cr, operator acting on the jth probe qubit and is the z component of the "total 
angular momentum" for all the qubits. 

Evolution under this Hamiltonian for a time t introduces a relative phase <p = yt between the two compo- 
nents of the superposition, changing the state of the probe qubits to {e~"^^^\0) + e"^^^\\))/ V2. The last set of 
Hadamard gates changes the parameter-dependent phases in the superpostion into amplitude (population) 
information. Thus the state of the probe qubits just before the readout is cos(^/2)|0) -i- sin(^/2)| 1 ). The final 
readout in Ramsey interferometry is done by measuring each of the qubits along the z-direction. This leads 
to a measured signal 
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FIG. 2: Quantum circuit for cat-state interferometry. 



The variance in the signal is 



{^^J,) = ^A^<aV,> = ^A^sin^^ 



(2.3) 



The uncertainty in the estimate of y from the measured signal in Eq. (12.2b is 

(aV,)!/^ 1 

(^TRamsey ^ ,,,,,,, i ^ ;= ■ (2.4) 

For the Ramsey Hamiltonian (12.11) . the generator of translations in y is tJ^. Thus, according to the 
quantum Cramer-Rao bound (11.131 ). the measurement uncertainty is bounded from below by 

1 1 1 



dj > 



(2.5) 



tm\ tN\\cTj2\\ tN 

The Ramsey interferometer described here does not achieve the best measurement uncertainty given by the 
quantum Cramer-Rao bound. The Hamiltonian //Ramsey that governs the evolution of the probe qubits is 
fixed by the choice of physical systems that are the qubits. Given a choice of probe qubits, however, we still 
have the freedom to choose an optimal initial state for the probe and an optimal measurement of the qubits 
to minimize the measurement uncertainty. It turns out that the best possible scaling for the measurement 
uncertainty can be achieved if the probe is initiahzed in an entangled, "Schrodinger-cat" state lioi 11]. 



B. Cat-state interferometry 

The quantum circuit that uses a probe initialized in a Schrodinger-cat state is depicted in Fig. |2l The 
Hadamard gate on the first qubit, followed by the controUed-NOT gates to the remaining qubits, initializes 
the probe in the state |cat> = (|0 . . . 0) -i- 11 . . . 1))/ V2. This state is often referred to as the Schrodinger-cat 
state because when the number of qubits is large, it is a superposition of two macroscopically distinct states. 

The probe qubits evolve under the same Hamiltonian (12.11 ) as in Ramsey interferometry. The 
parameter-dependent evolution of the probe for a duration t changes the probe state to (e~'^'''^^|0 . . .0) -i- 
^iNipii^Y ... 1»/ V2, where ip = yt. After the parameter-dependent evolution, one way to handle the readout, 
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depicted in the circuit above, is to subject the qubits to a sequence of gates that kick the phases picked up 
by the two components of the cat state into amplitudes on the first qubit, so that the state of the probe just 
before readout is [cos(A^^/2) |0> + sin(A^^ /2) 1 1 >] 10>®(^" ^ \ The readout of the probe can then be performed 
by measuring the cr^ operator on the first qubit. This leads to a measured signal and variance given by 

(o",;i> = cosA/^^ and {^cr^-^i) = ?,vc^ Nf . (2.6) 

The frequency of the y-dependent fringe in cat state interferometry is times greater than the frequency of 
the signal in ordinary Ramsey interferometry. This leads to an enhanced sensitivity in the estimate of y in 
cat-state interferometry, which achieves the Cramer-Rao bound: 

<AV,.i>'^2 1 

' ' m ■ 

C. Heisenberg-limited metrology with linear Hamiltonians 

We can put our interferometry circuits in a general setting by considering the case in which the probe 
units are arbitrary systems and the probe Hamiltonian is of the form 

N 

H]meai = y/llinear = y'^hj . (2.8) 

7=1 

Here the operators hj denote identical couplings to the probe units; the use of independent couplings to the 
parameter is the source of our appelation "linear" for this Hamiltonian. The generator of translations in y is 
K{y, t) - f/ziinear, SO the quantum Cramer-Rao bound (11.131 ) on the uncertainty in a determination of y takes 
the form 

1 1 

5y > ^ , (2.9) 

^ ~ fll^nearll tNiA - A) ' 

where A and A are the largest and smallest eigenvalues of the single-unit operators hj. Achieving the 
Cramer-Rao bound only requires using two levels of each unit, the eigenstates |A) and \A} corresponding to 
the largest and smallest eigenvalues of the operators hj, so we can always regard the units as qubits with 
|0>-|A>and|l>-|^>. 

The quantum circuit that represents a measurement protocol of this sort is drawn in Fig. [3] The dashed 
boxes highlight the three stages of this protocol: probe preparation, dynamics, and readout. All the probe 
units begin in a standard state |5 ). The arbitrary unitary operator P can then prepare any initial state as 
input to the dynamics. In the dynamics stage, the gates U^p imprint information about the parameter on the 
probe. The final readout stage includes an arbitrary unitary interaction R among the probe units and with 
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FIG. 3: Quantum circuit for a general linear interferometer 



an arbitrary ancilla system. This unitary followed by measurements on each subsystem in a standard basis 
can be used to perform any quantum measurement. The quantum Cramer-Rao bound (12.91 1 applies to all 
circuits of the above form. Indeed, the bound actually applies to somewhat more general situations in which 
the unitary operator R is interleaved with the gate dynamics and the results of ancilla measurements are fed 
back onto the probe 171]. 

If the preparation unitary P is omitted from the circuit, making the input to the dynamics a product state, 
then the uncertainty in the generator of y displacements is bounded by (A^K)^^^ < f V^(A - A)/2. The 
resulting bound on measurement uncertainty, from Eq. (11.121) . is 

1 



Sy > 



= <5yQNL ■ 



(2.10) 



t VyV(A - A) 

This bound, a general form of that for standard Ramsey interferometry, is called the quantum noise limit 
(QNL) or the shot-noise limit. The optimal 1 /t sensitivity for product-state inputs can be achieved 
by using initial state \S) = (|A) -i- \A))/ V2 for each unit and by making a final product measurement of 
an equatorial-plane spin component on each unit (in the qubit Bloch sphere formed from |0> = |A) and 
|1> = \A)). 

One can achieve the Cramer-Rao bound ( 12.91 ) by operating the circuit in a way that takes advantage of 
entangled input states. The preparation operator is chosen to take the initial product of standard states to the 
"cat-like" state (|A, . . . , A) + \A, . . . , A))/ V2. In the dynamics stage, this "cat-like" initial state is subject to 
a period of parameter-dependent evolution that changes it to (e~'^^'''|A, . . . , A) -i- e''^'^'^\A, . . . , A))/ V2. The 
readout process kicks back the differential phase shift into amplitude information, which produces fringes 
with frequency proportional to A'^(A - A), thus achieving the optimal measurement uncertainty, 

1 



6y = 



HL 



(2.11) 



tNiA - A) 

of the Cramer-Rao bound (12.91 ). This optimal measurement uncertainty, a general form of that for cat-state 
interferometry, is often called the Heisenberg limit. 
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FIG. 4: Quantum circuit for a general nonlinear interferometer. 



The general quantum-metrology scheme considered in this subsection indicates that probe preparation 
gives an enhancement of 1 / over the case where the probe qubits are initialized in a product state. 
Readout has already been optimized to take advantage of this entangled input, so we conclude that when 
the parameter-dependent dynamics acts independently on the probe qubits, Heisenberg-limited scaling is 
indeed the \IN scaling. The one remaining way o f exp loring whether the \/N scaling can be improved is 



to consider more general dynamics ||7 
subsection. 



mmg way or exptoring \ 



1711 : we turn to that possibility in the next 



D. Heisenberg-limited metrology with nonlinear Hamiltonians 

A generalized quantum-metrology scheme in which the dynamics of the probe is generated by a Hamil- 
tonian that includes all ^-body couplings between the probe qubits was first considered in [T]]. This nonlinear 
coupling Hamiltonian has the form 

^^nonlinear = y/jnonlinear = 7 hjj =y ^jl " ' " ^jk ■ (2-12) 

7=1 7'l.-.A=l 

The generator of translations in y is K(y, t) = f/inonimear> so the quantum Cramer-Rao bound for this dynam- 
ics is 

6y > —, , (2.13) 

fA^HAmax - Amin) 

where A^ax and A^in are functions of A and A, the largest and smallest eigenvalues, respectively, of the 
single-unit operators hj. For instance, if both A and A are positive, then A^ax - A*^ and Amin = A^ for all 
values of k. The other possible signs of A and A are discussed in [ 8] ; they all lead to a 1 /N'^ scaling. 

The quantum circuit for metrology with nonlinear Hamiltonians has the form shown in Fig. |4l This 
circuit has the same overall form as that for linear quantum metrology, with the same three stages high- 
lighted by the dashed boxes. The only difference comes in the dynamics stage, where the gate that imprints 
information about the parameter on the probe involves simultaneous couphng to all the probe units. 
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To achieve the 1 ///^ scaling made available by using a nonlinear Hamiltonian, the probe units have to 
be initiahzed in an entangled state that is very much like a cat state. Experimental limitations up till now 
haveprecluded making such cat-like states for large numbers of systems. To avoid this difficulty, Boixo et 
al. 180 analyzed the performance of quantum-metrology protocols employing nonlinear Hamiltonians when 
the initial state of the probe is a product state. In this case the optimal measurement uncertainty scales as 

The factors multiplying this scaling depend on the particular nonlinear coupling Hamiltonian |l8|]. It is 
noteworthy that the optimal 1 /tN'^^^^^ sensitivity can be achieved using product measurements of equatorial 
spin components in the effective qubit space formed from |A) and \A}. The key point is that for a /c-body 
coupling Hamiltonian, the use of a product-state input costs only a factor of ~ Va^ relative to the optimal 
sensitivity (12.131) . The quantum noise limit and the Heisenberg limit of linear metrology are a special case 
of this VA^ loss of sensitivity when using input product states as opposed to an optimal entangled state. 

With two-body couplings and an initial product state for the probe, a measurement uncertainty scaling 
as 1 /N^^'^ is possible. Since two-body coupling between all probe units is not an especially onerous re- 
quirement for a probe system, the prospect of improving upon the 1 /N Heisenberg scahng motivates us to 
investigate candidate systems for such metrology schemes. In the next section we consider a Bose-Einstein 
condensate (EEC) as such a candidate system with the aim of developing a detailed, realistic, and viable pro- 
posal for an experiment that achieves better than 1 /N scaling for the measurement uncertainty in quantum 
single-parameter estimation. 



E. Role of entanglement 



The more general point of view provided by nonlinear quantum metrology allows us to see exactly what 
benefit entanglement bestows on quantum metrology. Entanglement permits one to marshal the available 
resources in a quantum-metrology protocol into an initial state that can achieve the best possible scaling for 
the measurement uncertainty as laid out by the quantum Cramer-Rao bound. The use of an appropriately 
entangled input state purchases a sensitivity boost by a factor of ~ 1 / Va^ relative to the use of an optimal 
initial product state. Initial entanglement is, however, not necessary for getting to or improving upon the 
l/N Heisenberg scaling. 

When specialized to qubits, with hj = cr^- j/2, and to quadratic couplings, the nonlinear probe Hamilto- 
nian of Eq. (12.121 ) becomes //nonlinear - yJ?- This Hamiltonian generates entanglement during the dynamics 
stage of the protocol. Despite the evidence from the quantum Cramer-Rao bound that entanglement only 
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helps in the initial state, one might reasonably ask whether this dynamically generated entanglement plays 
a role in improving upon the IjN scahng. The j} probe Hamiltonian was analyzed in detail in Isl]. The 
optimal initial product state is [cos(;r/8)|0> + sin(;r/8)|l>]®". Evolution under the j} Hamiltonian for a short 
time t «; followed by a measurement of Jy, gives a measurement precision 6y = l/N^^'^, which is the 
optimal precision for this Hamiltonian and initial state. As a consequence of the evolution, however, 
the probe qubits become entangled and suffer from an associated "phase dispersion" that makes the mea- 
surement uncertainty large for separable measurements when yt becomes large. Far from being an aid, the 
generated entanglement seems only to make it impossible to achieve the 1 /N^^-^ sensitivity using product 
measurements. 



In ilSll it was pointed out that in addition to the Hamiltonian, //noniineai- = y^J^, can also be used to 
get an optimal measurement uncertainty 1 /N^^'^ when the probe qubits all start off in the state (|0) + 1 1 »/ V2. 
The NJ^ Hamiltonian does not produce entanglement between the probe qubits, nor does it produce phase 
dispersion. Indeed, the NJ^ Hamiltonian acts like a linear coupling whose strength is enhanced by a factor 
of N. In this case it is clearly the dynamics alone that leads to the enhanced scaling for the measurement 
uncertainty, since there is no entanglement between the probe qubits at any stage in the metrology protocol. 

On physical grounds, the NJ^ Hamiltonian cannot be a fundamentally linear coupling whose strength is 
enhanced by addition of qubits, but rather must come naturally from quadratic couplings to the parameter. 



Such a coupling does appear in the Hamiltonian for a two-mode BEC, as was pointed out in jlSh . We 
introduce this BEC implementation in the next section and discuss it in some detail. 

III. BOSE-EINSTEIN CONDENSATE AS A QUANTUM PROBE 
A. Nonlinear BEC interferometry 

The many-body Hamiltonian for a dilute Bose gas consisting of atoms of mass m in a trapping potential 
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Jf 1 ^ 

■ VVfr -I- y(r)\i/'''\j/ + -g\j/'''i$r'''\jlr\jlr I , (3.1) 



V{r) at zero temperature, in second-quantized notation, is given by 

where '^\r) and \i/(r) are creation and annihilation field operators that obey bosonic commutation relations, 

[Vi/(r), v+(r')] = 6^^\r - r'), [Vi/(r), Vi/(r')] = [^\r), v+(r')] - , (3.2) 

and the coupling constant g, for a dilute gas in which the inter-particle spacing is much larger than the 
scattering length, is related to the s-wave scattering length a by 

g = . (3.3) 

m 
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In a zero-temperature BEC, to a very good approximation, all the atoms are in the ground state i/fN(r), 
which is the A^-dependent ground-state solution (normaUzed to unity) of the time-independent Gross- 
Pitaevskii (GP) equation for a trapping potential V(r) and a scattering term with coefficient g: 

_|lv2 + Vir) + giN - mNlA = tiN^N . (3.4) 



2m 

Here /ia? is the chemical potential. At this level of approximation, the expansion of the field operator in 
terms of modal annihilation operators can be truncated to just one term, 

Mr) = ^]v(r)a , (3.5) 

where a annihilates a particle with wave function i/^Ar(r). The number operator a^a for this single mode can 

be treated as the c-number A'^ because the number of atoms is a constant. The Hamiltonian then reduces to 
the c-number mean-field energy for this single mode, H = EqN + jgrj^NiN - 1), where 

Eo = J^/r |^|V^^|2 + y(r)|^;vl') (3.6) 

is the single-particle kinetic plus trapping energy, and the quantity 

77iv = J dr\il,N{rf (3.7) 

is a measure of the inverse volume occupied by the ground-state wave function. The product gr]}^, which has 
units of energy, is a scattering strength normalized by this effective volume. The average number density in 
the atomic cloud is A^t/at. 

So far we have assumed that all the atoms in the BEC are in a single atomic state, but as mentioned 
earUer, we want these atoms to be two-level systems, or qubits, in order for them to serve as the probe units 
in the quantum-metrology protocols we are interested in. We therefore consider two-mode BECs in which 
the atoms can occupy one of two internal states, labeled |1) and |2). These two states are typically hyperfine 
levels of the atoms. In practice, the atoms are cooled to form the BEC while they are all in the same internal 
state, and then an external field is used to drive transitions between the two levels to achieve the desired 
coherent superposition of atomic population between the two levels. The effect we are looking for is the 
difference between the integrated nonUnear phase shifts experienced by the two levels, the difference being 
due to the different scattering interactions experienced by the two levels. This differential integrated phase 
shift is detected by driving a second transition between the levels, which transfers the phase information 
into the populations of the two levels. 

For an initial analysis of this scenario in this section, we make three simplifying assumptions: 
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1 . The external field that drives the transitions between the two states 1 1> and |2> acts only for a short time 
compared to the phase-shift dynamics that leads to the estimate of the parameter we are interested in. 
We therefore treat these transitions as effectively instantaneous and do not include the driving field 
in the Hamiltonian. 

2. The collisions between the atoms are elastic. Thus the only allowed scattering processes are |1)|1) — > 
|1)|1), |2>|2> |2>|2>, and |1>|2> |1)|2>, with scattering coefficients ^n, g22, and gn, where 
gafi = 4nti^aa/j/m - g^a, with Greek letters used to label the internal states. 

These first two assumptions imply that the many-body Hamiltonian takes the form 

^ = Z /^^ • + + \YjSap^dr v;¥l¥„¥^ , (3.8) 

where Vj/Q,(r) is the field annihilation operator for internal state a. In writing this Hamiltonian, we assume 
that any energy splitting between the two internal states has been removed by going to an interation picture. 
Our third assumption is by far the most problematical of the three. 

3. The two modes retain the same spatial wave function ipNir) as they evolve. Since the atoms that 
form the initial BEC are all in the state |1), in the mean-field approximation they all share the spatial 
wave function ifjt^{r), which is the A'^-dependent ground-state solution of the time-independent GP 
equation (13.41 ) with scattering coefficient g\\. Immediately after the nearly instantaneous action of the 
external field, the wave function for both internal states is t^Nir)- We further assume that the second 
internal state is chosen so that it sees the same trapping potential V{r). Even though the two internal 
states have identical initial wave functions and experience identical trapping potentials, their wave 
functions will gradually become different, because of the difference in their scattering lengths. What 
we are assuming now is that the integrated nonlinear phase shifts that we are interested in accumulate 
on a time scale that is shorter than the time scale for the two wave functions to differentiate spatially. 
Thus, for the present, we take the two wave functions to be identical. We return to the question of 
the time scale for differentiation of the two wave functions at the end of this section, in Sec. 1111 Fl 

Using the third assumption, we can write the field annihilation operators as 

"^air) ^ i//N(r)aa . (3.9) 

Since the total number of atoms is fixed, we can treat the total number operator. 



tC't ^ T ^ , T 



(3.10) 
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as a c-number N. We can then put the two-mode Hamiltonian in the form 

H ^ EqN + ^T]N ^ ga/jcij^ala^a^ = Ho + yi t]n{N - 1)/^ + JiJInJ^ ^ (3-11) 

where £0 and rj^ are as in Eqs. (13.61) and (I3.7I ). The operator is defined by 

= -(a\a^ -^2^2) ' (3-12) 

and we have also introduced a c-number energy, 

Ho = EoN + ^ l^teii + g22) + gxij rjN^^ - ^(gn + g22)riNN , (3.13) 

which includes the common-mode part of the mean-field scattering energy. Finally, we define two coupling 
constants that characterize the interaction of the two modes, 

n = ^teii - 522) and 72 = ^(gii + 822) - 812 ■ (3.14) 

The Hamiltonian (13.1 II ) is often called the Josephson approximation. 

The common-mode energy Hq in Eq. (13.111) can be ignored because its only effect is to introduce an 
overall phase in the evolved state of the probe. In the other two terms, we have (A^ - 1)/^ and /? couplings, 
suggesting that we might be able to measure the coupling constants yi and 72 with an accuracy that scales 
as l/N^^^ with the number of atoms in the BEC. 

To see how this works out, suppose the first optical pulse puts each atom in a superposition ci|l> -1- C2|2), 
where ci and C2 can be assumed to be real (i.e., the first optical pulse performs a rotation about the y 
axis of the Bloch sphere). For short times, we can make a linear approximation to in the Josephson 
Hamiltonian; i.e., we can set = ((/,) -1- AJ^f - (/j^ + 2{J.)AJ^, with (/j = N{c\ - c\)l2. The hnear 
approximation amounts to neglecting the phase dispersion and corresponding entanglement produced by 
the term. We need not make any such short-time approximation for the {N - 1)7^ term. Up to irrelevant 
phases, the resulting evolution is a rotation of each atom's state about the z axis of the Bloch sphere with 
angular velocity 

Y[(A^ - l)ri + N{c] - cl)y2\ - ^^^^^^[yi + {c\ - c\)y2\ = , (3.15) 

where in the second form, we approximate N as N- 1. Under these circumstances, the BEC acts like a linear 
Ramsey interferometer whose rotation rate is enhanced by a factor of {N - Vfrji^ , leading to a sensitivity that 
scales as 1/ VA''(A'^ - l)r}N - 1 /N^^^Vn- If 72 = 0, the optimal initial state has ci = C2 = 1/ V2, but if yi = 0, 
the optimal choice is ci = cos(7r/8) and C2 = sin(;r/8) [Sl]. 
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Achieving a l/N^^'^ scaling requires that ijn have no dependence on A'^. As noted above, however, 77^^ is 
a measure of the volume occupied by the ground-state wave function 1^^. As atoms are added to a BEC, the 
wave function spreads because of the repulsive scattering of the atoms, thereby reducing tjn as A'^ increases. 
To pin down how the measurement accuracy scales with N, we need to determine how tjn behaves as a 
function of N. 

B. Two critical atom numbers 

Since we first create a BEC of N atoms all in hyperfine state |1), before putting them in a superposition 
of states |1) and |2), we can focus on the A/^-dependence of rj^ for a single-mode BEC of atoms in state 
Thus, in this subsection and the next two, we deal with the single-mode GP equation (13.41 ) with g = gn and 
a = ail. 

An obvious strategy to suppress the A^-dependence of rj^ is to constrain the BEC within a hard-walled 
trap so that it cannot expand as more atoms are added. BECs effectively confined to two or one dimensions 
and held in power-law trapping potentials along these dimensions are the sort found in real experiments. 
Thus we look at the dependence of t/at on A'^ for a BEC that is loosely trapped in d dimensions, referred to 
as longitudinal (L) dimensions, and tightly trapped in D = 3 - d dimensions, referred to as transverse (T) 
dimensions. We assume that in the longitudinal dimensions, the atoms are trapped in a power-law potential 
of the form 

VL{r) = ^kr'^, q=\,2..., (3.16) 

and that in the transverse dimensions, the trapping potential is harmonic, 

Vrip) = . (3.17) 

The parameter q characterizes the hardness of the longitudinal trapping potential. We deal with a 3D trap 
by setting D = 0, meaning there are no transverse dimensions. 

When A^ is small, the mean-field scattering energy is negligible compared to the atomic kinetic energy of 
the atoms and the trapping potential energy. In this situation, the scattering term in the GP equation can be 
neglected, and the ground-state wave function is the solution of the Schrodinger equation for the trapping 
potential Vi(r) + Vjip)- As more atoms are added to the BEC, the repulsive scattering term in Eq. ( 13.41 ) 
comes into play and causes the wave function to spread. We define two critical atom numbers, A^^, and A^j-, 
which characterize the onset of spreading in the longitudinal and transverse dimensions. The lower critical 
atom number, A^^, is defined as the atom number at which the scattering term in the GP equation is as large 
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as the longitudinal kinetic-energy term and thus characterizes when the wave function begins to spread in 
the longitudinal dimensions. The upper critical atom number, Nt, is defined as the atom number at which 
the scattering term is as large as the transverse kinetic energy and thus characterizes when the wave function 
begins also to spread in the transverse dimensions. The notion of an upper critical atom number only makes 
sense for ID and 2D traps and not for d = 3. 

For small atom number, i.e., A'^ «; Nl, as just noted, the scattering term in the GP equation can be 
neglected, and the ground-state solution of the GP equation is the A'^-independent, product ground state of 
the Schrodinger equation: 

Mp,r)^Xoip)Mr). (3.18) 
Herexoip) is the Gaussian ground state for the transverse dimensions, 

whose corresponding probability density has half-width 

\2ma>T/ 

and (poir) is the bare ground state for the loosely confined longitudinal dimensions. We can estimate the 
half- width of (po by equating the trapping potential energy and the kinetic energy (KE) per dimension, i.e., 
= fp-jlmr^, which gives 

. (3.21) 

In accordance with our assumptions, we assume that ro is much larger than po- A hard-walled trap in the 
longitudinal dimensions corresponds to the hmit q ^ oo with ro held constant. 

The trapped BECs we consider are thus characterized by three length scales: (i) the scattering length a; 
(ii) the bare transverse trap half- width po; and (iii) the bare longitudinal trap half- width tq. Typical values, 
which we use for estimates in the following, are a = 10 nm, po = 1 yum, and ro = 100 /im. For ^^Rb atoms 
(which have a - an - 5.3 nm), the corresponding transverse trap frequency is vt - 58 Hz; we can also 
identify an approximate longitudinal trap frequency, 

-- = ^-f A-lO-^Hz, (3.22) 

associated with the bare longitudinal ground state. 
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Whenever the wave function is a product of transverse and longitudinal wave functions, tjn is also a 
product, r]i^ - tjjtjl- When N «; Nl, i]n = '7o is independent of N since 

riT = f d''pkoip)f= ^l, n ' (3-23) 



VL = J A|0o(r)|'-T7^, (3.24) 



1 

where is the volume of a unit sphere in d dimensions (Vi = 2, V2 = n, and V3 = 4n/3), The lower critical 

atom number, A^^, is defined by setting 

if tf I a 

(longitudinal KE) - (scattering term) - {Nl - i)gm - Ti—i^L - 1) 77-7 , (3.25) 



where 

is a geometric factor (fii = \,^2 = V^/4,jS3 - 1/6). Thus we define 

NL-l^/3d-l-t . (3.27) 
a \ ro / 

For the typical length scales mentioned above, the lower critical atom number is about 1 700 for a 3D trap, 
45 for a 2D trap, and 2 for a ID trap. The small value of A^^, for a ID trap is the reason we retain the -1 
wherever it appears in our discussion of atom numbers, even though it could be dropped in most situations. 

For Nl < N ^ Nj, the tight confinement in the transverse dimensions means that the wave function 
continues to be a product, 

if/Nip, r) = Xoip)(f>Nir) , (3.28) 

but with the longitudinal wave function satisfying a GP equation, 

-^^1 + VLir) + g{N - l)?7rl0ivM = MicpN , (3.29) 

where Hl = IJ-n - Dfiojjil is the longitudinal part of the chemical potential. As atoms are added to the trap 
in this intermediate regime, the wave function spreads in the longitudinal dimensions. We can estimate the 
longitudinal half-width rj^ by noticing that r]j^ - tjttjl, where rjj is given by Eq. (13.231 ) and 

1 

and then equating the attractive longitudinal trapping potential energy (PE) to the repulsive scattering term: 
1 



TiL = Jdr \cl>N{r)f - , (3.30) 



-krl, ~ (longitudinal PE) 



Jf 1 a fp- N 



- (scattering term) - (A^ - \)gTiN - :^(A^ - 1)^-5-7 = ' ^^'^^^ 

2m Pd p^fj^ 2mNL- I 
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where we have used Eq. (13.271) in the last step. This leads us to define 



rN i N-l 



l/{d+q) 



ro \Nl - 1 

We now define the upper critical atom number by setting 



(3.32) 



-2 



Nt-i 4~ 

(transverse KE) - (scattermg term) - {Nj - ^)griNT = ' (3.33) 



where rj is the longitudinal half-width at the upper critical atom number, 



-\Ar ■ (3.34) 

'•o \A^L-1/ 

Using Eq. (13.331 ) and the definitions in Eqs. (13.341) and (13.271) we end up with the definition 

Nt-1^{Nl-\)\-\ =Pd-\-] . (3.35) 

\Po / a\pol 

We stress that the notion of an upper critical atom number only makes sense for ID and 2D traps and not 
for d = 3. Using the typical values mentioned above, we have that the upper critical atom number for a 
harmonic longitudinal trap (q - 2) is about 4 x 10^ for a 2D trap and about 10^ for a ID trap; for a hard 
longitudinal trap (q oo), Nt is about 4 x 10^ for a 2D trap and about 10^ for a ID trap. Using Eq. (13.351) 
we can rewrite the longitudinal radius in Eq. (13.321) as 



2/? / 1 \y(d+q) 



(3.36) 



ro \po/ \A^7- - 1 
It should be noted that 

= ^ • (3-3^) 

Po \Po Pd I 

For a ID trap, this gives rj - a{NT - 1), making the relation between rj and Nj independent of the 
parameters of the trap. Another way of thinking about Eq. ( 13.371 ) is that the number density at the upper 
critical atom number, 

A^r 1 

TTl " ' (3.38) 

is independent of the properties of the longitudinal trap, with typical value 10'^ cm"^. 

As the atom number increases from Nt, the transverse kinetic energy becomes unimportant compared 
to the transverse trapping energy and the scattering term. The wave function continues to spread in the 
longitudinal dimensions and also spreads in the transverse dimensions, with the longitudinal and transverse 
radii, rj^ and pi^, given by 



^^r^ - '^'^^tPn - (scattering term) 



(A^ - DgrjN ^{N- 1)^— -1— , (3.39) 

Vop'fiVdri, 
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which leads us to define in the regime A'^ » Nt, 



^^(^^) , (3.40) 
ro \2poPo/ 



C. Renormalization of the nonlinear interaction terms and the sensitivity scaUng 

The estimates in the previous subsection tell us how tjn scales with atom number. For atom numbers 
smaller than the lower critical atom number, tjn has the constant value ijq, a consequence of the fact that 
the repulsive scattering has negligible effect on the atomic wave function. In the intermediate regime of 
atom numbers, i.e., for atom numbers between Nl and Nt, the wave function expands in the longitudinal 
dimensions, making rjj^ scale as 

1 /M ]\dl(d+q) 

For atom numbers above the upper critical atom number, as the wave function spreads in all dimensions, 77 a? 
scales as 

In the measurement schemes we contemplate, the uncertainties in determining y\ and 72 scale as 

(5ri.2 ~ ^ ^ , (3.44) 

^^N{N-\)7^N 

where in the final form we neglect 1 compared to N. For atom numbers below A^^, the scaling exponent ^ is 
3/2; for A^L «: A/^ «; A^r, it takes on the value 

lA'^IlISL; (3.45) 

2 a + q 2{a + q) 



and for A'^ » Nt, ^ is given by 



3_3-d + 2d/q 
^ 2 5-d + 2d/q 



For atom numbers above A^7, harmonic ID and 2D traps have ^ = 9/10, a hard-walled ID trap has ^ = 1, 
and a hard-walled 2D trap has ^ = 7/6. Our main interest is the intermediate regime of Eq. ( 13.451) . The 
scaling exponent in this regime is plotted in Fig.|5]as a function of q for ID, 2D, and 3D traps. 
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FIG. 5: (Color online) Sensitivity scaling exponent^ = (d + 3q)/2(d + q) in the intermediate regime of atom numbers, 
Nl N Nt, plotted as a function of hardness parameter q of the longitudinal trapping potential for ID (red, 
solid), 2D (green, dashed), and 3D (blue, dotted) traps in the intermediate regime of atom numbers (for 3D traps, 
there is no upper critical atom number A^j-). To achieve super-l/A^ scaling (^ > 1) requires q > d. A harmonic ID 
trap has super- scaling ^ = 7/6, a harmonic 2D trap has Heisenberg scaling ^ = 1, and a 3D harmonic trap has 
sub-l/A^ scaling ^ - 9/10. This sub-l/A^ scaling for 3D harmonic traps is still markedly better than the QNL scaling 
of ^ = 1/2. A hard-walled trap {q — > oo) in any dimension has ^ = 3/2. 

D. Thomas-Fermi approximations 



Although we have determined how the scaling exponent behaves with d and q, we can do a better job of 
evaluating rjf^, determining more precisely the constants in front of the scaling, by using the Thomas-Fermi 
(TF) approximation. In the intermediate regime of atom numbers, the wave function is the product (13.281 ). 
with the longitudinal wave function (/>iv(r) satisfying the GP equation (13.291 ) in d dimensions. When N is 
much larger than Nl, we can ignore the kinetic-energy term in the reduced GP equation, which gives the TF 
probability distribution. 



2 f^L-krVl 
f>N(r)\ = 



(3.47) 



(A^- l)g7]T 

Since \(pi\i(r)\^ must be positive, the radial extent of the BEC in the longitudinal dimensions is bounded by 
such that 



Normalization yields 



-f 



d^rlMrr = hiN,d,q) 



(3.48) 



(3.49) 
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Integrals over TF probability distributions, such as /i , are defined and evaluated in an Appendix. Using 
Eq. (lAT/l ). we find 



r/v ( d + q N - I 



[cf. Eq. (13.321 )1 or, equivalently, using Eqs. (lA.ll ) and (IA.4I ). 



l/{d+q) 



lUL d + q I 1 (d + qY"'^''^ 



(3.50) 



(3.51) 



(A^ - l)gr]T q Vdrfj Vd^X q ) \N-l 
Now, from Eq. (IA.8I ). we can find 

''^ = ^^^^'^'^^ ^ ..2,(A.-1),,, ^ ^^l— ) I— ) ^^-^'^ 
and thus determine tjn = tjttjl [cf. Eq. ( 13.421 )1. Numerical computation of rj;^ in the intermediate regime 
indicates that this expression is quite accurate in spite of the approximations that went into obtaining it. 

When N is much larger than Nt, we can again use a TF approximation, this time ignoring the kinetic- 
energy term in the 3D GP equation (13.41) . which gives the probability distribution 

T Um - mcol.p^/2 - kr^ 12 

\^N{p, rf . (3.53) 

(A^ - 1)^ 



Positivity of this distribution requires that p < pn and r < fj^ip), where 

-mo?jp]^ = fiN , (3.54) 



2 



^kf^ip) ^ ^ma»|(p^-p2). (3.55) 

The extent of the atomic cloud in the longitudinal direction is characterized by r^? = ri\i{0), i.e., = 
[(ro/2po)(PAr/po)]^^^- Normahzation yields 

1 - Jd^p d'^r \il/N{p, r)|2 = ^i(A^, d, q) . (3.56) 

Using Eq. (lA.UI) . we find 

(pnX''''^^'^^'' 1 4(47r)^/222^/, ^ _ 1 



ypoj dJi+diq{D,2)J\{d,q) Sd-i Nt - I 

[cf. Eq. (13.411 )1. We can also write 

UN ^ 1 (Pof" (pA' Nt-1 

{N-l)g 4i4n)D/2pDvy^\rol [poj N-l 



(3.57) 



(3.58) 



and 



The scahng of t/at agrees with that in Eq. (13.431 ). 
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E. Bose-condensed '^^Rb atoms 



A good candidate for implementing the generalized metrology protocol is a BEC made of Rubidium 



(^^Rb) atoms 1231 



most experiments 



24i|25L|2^. Initially the atoms in the condensate are in the ground electronic state. In 
25L|2a], the = 1; Mp = -1> = |1) state is trapped and cooled to the condensation 
point. Once the atoms in |1) have accumulated in the condensate ground state, a two-photon drive is used 
to couple the |1) state to the \F = 2; Mp = +1) = |2> state. The two-photon drive involves applying 
both microwave and radio-frequency electromagnetic fields to the condensate. Both hyperfine states are 
not cooled simultaneously to form a condensate in a superposition of the two states because the lifetime 
of atoms in the |2) state in a trap is much shorter than the lifetime of atoms in the |1) state. As mentioned 
earlier, we assume that the driving field that initializes the atoms in a desired superposition of the 1 1 ) and |2) 
states is instantaneous in comparison to the dynamics that is part of the parameter-estimation process. 

The s-wave scattering lengths for the three processes, |1>|1> |1>|1>, |2>|2) |2>|2>, |1>|2> |1>J2) 
are nearly degenerate for ^^Rb, with the ratios {022 '■ «i2 '■ <^ii] - {0.97 : 1 : 1.03) (an - 5.31 nm) 12311 . 
These values for the scattering lengths mean that 72 = (gii + 822)/'^ - gn - for *^^Rb. Therefore, the 
^^Rb BEC realizes the generalized quantum-metrology protocol with just the 'Xi(A'^ - 1)??^/^ coupling in the 
Hamiltonian (13.111 ). 

The optimal initial state for metrology with the yiiN - l)!]^!^ coupling is the one in which all atoms 
are initialized in the equatorial plane of the Bloch sphere, say, in the -i-l eigenstate of cr^, (|1> -1- |2))/ V2. 
The quantity estimated by this measurement scheme is yi = (gn - 522)/2, which is small, but finite in the 
case of a ^^Rb BEC. Once the probe is initialized in the optimal initial state, we let it evolve for a time t 
under the influence of the yiiN - \)r]j^J^ Hamiltonian, which simply rotates the state of each atom about the 
z axis of the Bloch sphere with angular velocity yiiN - \)r]j^lfi. At the end of this evolution, we measure 
an equatorial component of J {Jy for short evolution times), which is achieved by a njl pulse about the 
desired equatorial axis, followed by a measurement of 7^, i.e., of the difference in the populations of the two 
internal states. 

Precision experiments with two-component ^^Rb BECs in modestly nonspherical, harmonic traps have 
been reported in l27i and |28]. These experiments were carried out with atom numbers in excess of 100000 



and thus operated in the full TF regime well above the upper critical atom number. 
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F. Differentiation of the spatial wave functions 

The strongest assumption we made in obtaining the Hamiltonian (13.111 ) was that the wave functions for 
the two modes remain identical throughout the duration of the proposed measurement scheme. Here we 
examine this assumption more carefully. 

We noted earlier that even if the two modes see the same trapping potential, the difference in their 
scattering lengths will cause the two wave functions to evolve differently [24]. The initial effect of the 
difference in scattering lengths is to produce a relative phase between 1 1) and |2). This relative phase depends 
on the local density within the condensate. The integrated (or average) part of the relative phase provides 
the signal for our measurement protocol, whereas the residual position-dependent part of the relative phase 
reduces the visibility of the fringes on which the signal relies. For our protocol to succeed, we need the 
integrated phase to accumulate more rapidly than the residual position-dependent phase. Yet a further effect 
is that the position-dependent phases drive differences between the atomic densities associated with the two 
hyperfine levels, but as this occurs on a longer time scale than the accumulation of the position-dependent 
phase shift, we do not consider it here. 

We can analyze this scenario in the following way. Initially all atoms are in the state i/'/v(p, r)|l). After 
the first optical pulse, the state becomes if/Nip, '")(ci|l) -i- C2|2)), where ci and C2 ai^e the amplitudes to be 
in the hyperfine states. We can assume that ci and C2 are real, i.e., that the initial optical pulse produces 
a rotation about the y axis of the Bloch sphere. The different scattering lengths make the wave functions 
for the two modes evolve differently, so that after a time t, the atomic state becomes ci(Aiv,i(p, r, OU) + 
C2^N,2^p, f, t)\T), where the wave functions for the two modes evolve according to time-dependent, coupled 
GP equations: 

^iV,a . (3.60) 



in- 



dt 



V2 + V + (^ - 1) 2 



The second optical pulse is a n/2 pulse about an equatorial axis of the Bloch sphere. For the discussion 
here, we assume that this rotation is about the x axis so that subsequent counting of the populations of the 
two hyperfine levels is equivalent to measuring Jy before the second optical pulse. The state after the pulse 
is 

111 / 

Cil/^ATj— (|l>-/|2»-l-C20'iV,2— p(-2ll> + |2» = —{cnl/N,l-iC2f/^N,2)\^} f(ci iAa',1 + 'C2lAiV,2)|2) . (3.61) 

V2 V2 V2 V2 

The corresponding probabilities to be in the two states, 

Pia = + 2ciC2lm((l/rA,_2|lAiV,l»] , (3.62) 
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are determined by the overlap of the two spatial wave functions, 



(3.63) 



In a ground-breaking set of experiments, Anderson et al. [28] measured the position-dependent phase 
shifts in a two-component ^^Rb BEC, trapped in a modestly nonspherical trap, and saw the associated re- 
duction in fringe visibility. The details of the experiment were shown to be well accounted for by numerical 
integrations of the two-component GP equations (13.601 ) with a loss term included. The experiment was 
carried out with atom number A'^ - 1.5 x 10^, well above the upper critical atom number. 

To compare the time scales for the integrated and position-dependent phase shifts in our protocol, we 
assume that we are operating in the intermediate regime of atom numbers, i.e., Ni < N <^ Nj- In this 
regime, the wave functions for the two modes factor into transverse and longitudinal wave functions, i.e.. 



where xo is the time-independent, Gaussian ground state in the transverse dimensions and the longitudinal 
wave functions obey time-dependent, coupled, longitudinal GP equations. 



To estimate the time scales, we assume that is large enough relative to the lower critical atom number 
to justify the TF approximation in the longitudinal dimensions, thus allowing us to ignore the kinetic-energy 
terms in the coupled GP equations. With these assumptions, the probability densities do not change with 
time, i.e.. 



(3.64) 




(3.65) 



(3.66) 



and the evolution under the coupled GP equations only introduces a phase. 




(3.67) 



This gives an overlap {4'n,2\4'n,\) - {(f>N,2\(f>N,i} = j d^^qoe "'^^^^ where the relative phase is given by 




(3.68) 



with 



= Ci(gii - gn) - cl{g22 - gu) = n+ (ci - cl)y2 ■ 



(3.69) 
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In the second equality of Eq. (I3.68I ). we have separated out the integrated phase shift, which has angular 
frequency 



[cf. Eqs. (13. isl and (13.221 )1. leaving the residual position-dependent phase shift as a correction. The final ex- 
pression for QjM uses the TF approximation to evaluate tjn in the intermediate regime. For the ^^Rb protocol 



fringe visibility in Eq. (13.621 ). 

It is worth emphasizing how this approach based on coupled GP equations differs from use of the Joseph- 
son Hamiltonian (13.111) . Although the GP equations yield a position-dependent phase, which cannot be 
obtained from the Josephson Hamiltonian, this comes at a price: the integrated relative phase in Eq. (13.681 ) 
amounts to making the linear approximation to /? described in the paragraph containing Eq. (13.151 ). The 
linear approximation is essential because the coupling does not preserve product states, whereas the GP 
equations assume a product state. It means that the GP equations miss the phase dispersion generated by 
the coupling and the associated dynamically generated entanglement. 

We can now write the overlap as 

(MM = e"'""' j d'rqoe-'''^'^"^"'' - ^>-'"-' exp / d'rqoiqo - r^tf'j , (3.71) 

where the second expression comes from expanding the exponential inside the integral to second order and 
then converting to an equivalent Gaussian at the same order. The contribution from the first-order term 
vanishes since r]L - j d'^ql- We can identify a time scale Tpd for the position-dependent phase as the time 
set by the half-width of the Gaussian, i.e.. 



The final form comes from using the TF approximation (13.471 ) for the density and the results in the 
Appendix to evaluate the integral. 

What this result means is that to retain good fringe visibility, our protocol will generally be restricted 
to operating well within the first fringe. One can expect, however, that as the longitudinal trap becomes 
more hard-walled, the TF density becomes more and more fiat-topped, eventually approaching a box, with 
the result that the residual position-dependent phase shift becomes smaller and smaller. This expectation is 
borne out by Eq. (13.721) . which reports that Tpd gets larger as the hardness parameter q increases; e.g., for a 
ID trap with q = 10, OAfTpd - 8. 




(3.70) 



outlined in Sec. IIII El in which 72 is essentially zero, we choose c? 



= C2 = 1 / 2 in order to maximize the 




(3.72) 
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To investigate further this way of reducing the effect of the position-dependent phase requires numerical 
simulations and more accurate approximation procedures, both of which we have undertaken. Initial results, 
to be reported elsewhere, suggest that things turn out better than is suggested by the crude approximations 
that go into Eq. (I3.72I ). 



G. Other practical considerations 

1. Loss of Atoms 

For a one-component BEC, the total number of atoms for a given trap is Hmited by three-body losses. 
This process is usually the most significant loss channel, with all other losses being negligible. For a 
two-component BEC, however, things are different because other loss channels, such as inelastic two-atom 
(spin-exchange) collisions, become significant even when the number of atoms in the trap is such that three- 
body collisions are unimportant. Just as in the case of three-body collisions, the spin-exchange collisions 
can be considered as a process that leads to loss of atoms from the trap. 

Spin-exchange collisions in a two-c omp onent BEC of ^^Rb atoms in the two hyperfine levels we are in- 
terested in were considered in ll27ll and 12811 . The effect of inelastic spin-exchange interactions was modeled 
by including non-Hermitian potentials in the coupled GP equations (13.601 ): 

~{N -\)Tncl\i]/2?- for model, (3.73) 
-|(A^ - l)(ri2C>i|2 + r22c2|^2l') for mode 2. (3.74) 

The loss constants in ^''Rb were measured to be ri2 = 0.780(19) x 10^'^cmVs and r22 = 1.194(19) X 
10~'^cm^/s. If we assume that the wavefunctions are the same for the two hyperfine states, as in the 
short-time analysis of Sec. IIII F[ the integrated effect of the spin-exchange losses across the atomic cloud is 
characterized by a decay constant 

^^ (jv-i);7^(ri2 + r22c^) ^^^^^ 

We can get an idea of the importance of spin-exchange losses by comparing F to the angular frequency O^f 
for the integrated phase shift. The ratio of interest for comparing coherent and decoherent processes is thus 

F /i(Fi2 + F22/2) m F12+F22/2 1 



Q.N 2yi AnTi a\i - a^i 19 



(3.76) 



where we specialize to the case c\ = c\ = Xjl relevant to the ^^Rb protocol and the final estimate applies to 
that protocol. This ratio indicates that the proposed protocol can obtain an estimate of y\ with better than 
1 I'M scaling before atom losses degrade the sensitivity. 
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It is also to be noted that an advantage of using a measurement scheme that uses product states is that 
loss of atoms from the BEC does not change the sensitivity scaling, since loss of particles from a product 
state does not damage any coherence. There is a decay in the signal strength given by a factor e~^', which 
would require us to complete the experiment before too many atoms are lost, but the ratio (13.761 ) provides 
a window for doing this. The discussion in Sec. IIII Fl suggests, however, that differentiation of the spatial 
wave functions for the two modes becomes a limiting factor on the duration of the experiment before loss 
of atoms becomes an important consideration. 

2. Number uncertainties 

In real experiments the number of atoms in a BEC is not known to arbitrary precision as we have assumed 
so far. Thus we have to consider what happens when the number of atoms in the BEC is not fixed from trial 
to trial. 

To analyze this situation, let p(No) denote the probability that the number of atoms participating in our 
measurement protocol is A^o- The final step in the protocol is to count the number of atoms in the two 
hyperfine levels. The difference between the two counts is used to estimate the parameter, here denoted as 
y; the sum can be used to refine the estimate of the number of atoms that participated in the protocol. 

We let N[ and N'2 be the number of atoms that would be counted by an ideal counting procedure. We 
generally work in terms of the total number of atoms, A^o - N[ + N'2, and the difference, m' = (N[ - N'^)!!, 
normalized by a factor of two to match the eigenvalues of J^. Quantum mechanics gives the conditional 
probability q{m'\NQ,y) for a measurement of 

The counting is not completely precise, so we introduce independent conditional probabilities, p{N\\N[) 
and p{N2\N'2), for counting A'^i and N2 atoms in the two levels, given the ideal counts. We can think of these 
two probabilities as describing processes in which condensate atoms are missed or non-condensate atoms 
are counted by mistake. In addition, in a complete analysis of the protocol, we would need to include the loss 
of atoms, discussed in the previous subsection, in this analysis. As already noted, we are mainly interested 
in the total number of atoms counted, N = N\ + N2, and the normalized difference, m = (Ni - A^2)/2. In 
the absence of a better model, we assume, to illustrate the effect of number uncertainties, that p{N\\N[) 
and p{N2\N'2) are independent Gaussian random processes, with mean Nj = Nj and variance A^A^y - cr^. 
Under this assumption, N and m become independent Gaussian random processes, described by conditional 
probabilities p{N\Nq) and p{m\m'), which have N = Nq, m = m', A'^N = Icr^, and A^m - cP' 12. 

The probability this model gives us directly is 

p{N, m, m',No\y) = p{N\No)p{m\m')q{m'\NQ, y)p{No) , (3.77) 
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The probability we need in order to evaluate the sensitivity of our protocol is the conditional probability for 
m, given the parameter y and the measured total number of atoms, A'^: 

p{N, m\y) 



p{m\N,y) = 



p{N\y) 

^ p(N,m,m',No\y) 



m' ,Nq 



p{N,m,m,No\y) 

m,m',No 

= P(fn\m')q{m'\No,y)p{No\N) . (3.78) 

m',No 

In the final form, p{No\N) = p(N\No)piNo)/ p{N) is the conditional probability for A^o atoms to have partic- 
ipated in the protocol, given the measured total count A'^. It quantifies the refinement in the knowledge of 
A'^o provided by the total count. 

The quantities that go into determining the sensitivity are the mean and second moment of m, calculated 
from the probability (13.781) . 

mN,y = Y,<^^,)No,yPiNQ\N) , (3.79) 



No 

{m^)N,y = \(t' + Yj{(Jz)lo,y + i^^ JdNo,y)pmN) , (3.80) 



where {J7}No,y = ("^O/Vo,/ ^^id (A^/j)Aro,y - (A^mOw^.y are the mean and variance of X calculated from 
the quantum-mechanical probabilities. If cr is much less than the initial uncertainty in A'^o, which is itself 
somewhat less than A'^o (depending on the care taken in loading the trap), then the measured total count 
A'^ gives a very good, improved estimate of the number of atoms that participated in the protocol; under 
these circumstances, the probability p{No\N) is peaked at the measured value A'^, with half-width given very 
nearly by a. The quantum-mechanical expectation values vary over a range from -No/2 to +No/2, so as 
long as o" «; A'^, we can evaluate the averages over p{Nq\N) at the mean value A'^ with little error, thus giving 
mean m^^y - {Jz)N,y and variance (A m)^^^ = cr /2 + (A J,)N^y. The resulting measurement uncertainty in 
determining y, 

2 ^ ^ ^ ^ -J£l± , (3.81) 

\dmN,yldy\^ \d{J,)j^,yldy\^ 

has the quantum-mechanical scaling and nearly the optimal sensitivity, provided we can count atoms to 
better than Va^, i.e., cr < V^. Ultimately, what this result expresses is that the variance of the measurement 
of in our protocol is of order V^V, so we need to know the number of atoms to this same accuracy. 
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IV. CONCLUSION 

This paper serves two purposes. The first is to extend the discussion of Heisenberg-limited quantum 
metrology from its traditional focus on a 1 /A/^ scaling for measurement uncertainty. Our discussion centers 
on the role of the dynamics of an A/^-qubit quantum probe in determining the quantum Cramer-Rao bound 
for single-parameter estimation. Looking at quantum metrology using the language of quantum circuits 
makes it easy to see that abandoning the usual independent couplings of the parameter to the qubits in favor 
of nonlinear couplings can yield scalings better than 1 IN. With ^-body couplings, it is possible to achieve 
sensitivities scaling as \IN^. Although the IIN'^ scaling requires entangled input states, 1/N''~^^^ scalings 
can be obtained with initial product states. Thus a sensitivity scaling as 1 /N^^^ can be achieved if quadratic 
couplings to the parameter can be engineered; moreover, particular quadratic couplings yield this sensitivity 
even though the state remains unentangled under the dynamics, thus showing that super- 1 /N scaling can be 
achieved without any entanglement. 

The second purpose of this paper is to show that a two-component BEC is a promising candidate system 
for a proof-of-principle experiment that demonstrates scaling better than 1 /N. A simplified analysis of the 
system, based on strong assumptions, but followed by a more detailed analysis of the realm of applicabil- 
ity of those assumptions, shows that such an experiment might indeed be realizable. This work motivates 
further, yet more detailed analyses and numerical simulations of the experiment. We have undertaken such 
further investigations of the proposed metrology scheme, and this further work, to be reported elsewhere, 
supports the conclusions reached in this paper. Our numerical studies include computing the ground-state 
solution of the time-independent GP equation for different values of A'^ in order to find the exact depen- 
dence of T]]\i on N. Numerical integration of the time-depenent, coupled, two-mode GP equations (13.601 ) 
is then used to compute the expected signal (13.621 ) in order to compare it with the theoretical prediction in 
Eq. (IXTTl) . 

The quantity that is measured in the proposed metrology protocol is essentially a constant. Estimating a 
constant using sophisticated quantum measurement schemes is interesting only as a proof of principle, be- 
cause there is nothing to preclude estimating the same constant using much simpler, classical measurement 
techniques. Since the measured quantity is a constant, we have the time to perform whatever number of 
repetitions of the simplest estimation procedure is required to achieve the desired accuracy. Metrology pro- 
tocols of the type described here are relevant in circumstances where there are constraints on the available 
time or on the available number of qubits. The available time can be constrained, for example, because the 
quantity that is being measured is changing, as in the case of gravitational-wave detection or magnetometry. 
There can be further time constraints placed by decoherence of the probe qubits. In such scenarios, picking 
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the optimal metrology scheme with the best measurement uncertainty, given the constraints, is of primary 



importance 12911 . For our proposal using a BEC, one possibility is to work around a broad Feshbach reso- 
nance that makes the scattering lengths sensitive to external magnetic fields. We might then be able to use 
our scheme for high-precision magnetometry. 
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APPENDIX: INTEGRALS OVER THOMAS -FERMI DISTRIBUTIONS 

In the intermediate TF regime, i.e., Nl ^ N ^ Nt, vj& need to do integrals over the TF probability 
density (13.471) . 

Il(N,d,q) = J d^'rl^Nirf 



k/2 


{N 






k/2 


{N 








(N 


- l)CTr 



/j^''^Sd-\_Ji{d, q) 



~4Sd-iJi{d,q) , (A.l) 



where Sd-\ = dVd is the area of a unit sphere ind - I dimensions and 



J,{d,q)^ rduu'-\\-u'^y = - . " rJ/'^^^ (A-2) 

Jo sm(Z;r)r(-/) qY{dlq + 



for/> -1. 

It is easy to see that 

Jx+i{d,q) {x+\)---{x + l) 



Jxid, q) (d/q + x + 1) • • • (d/q + x + I) 
Combined with Jq = I Id, this gives, when / is a nonnegative integer. 



(A.3) 



/! q^ 

" d{d + q){d + 2q)---{d + Iq) ' 
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Notice also that 



W, 2) = dv sin'^-V cos^'+i V = IfJ,^}^^^;' . (A.5) 
Jo 



2T{dl2 + 

Now we use 



Sd^i = , (A.6) 



(A^-l)CTr /^^^ N-l 



to write 



/i(A^,t/,^) = — — — - = — — — — — - (A.7) 

\ro/ d + q\roj N-l 

and 

I,{N, d, q) = h{N, d, q)^^Y^^ ("rp^^f^ ) ■ (A.8) 

Ji{d,q) \{N - Vfgmj 

In the upper TF regime, i.e., » Nt, we need to do integrals over the TF probability density (13.531 1. 
Ki{N,d,q) = J J^pA|(Aw(p,r-)P' 



k/2 



iN-l)g 



mojjl2 \ [ mojj ' 



{N-\)g} \ k 

I I..., ,2 



p'^'^^^^'^'^'''^SD-\Sd-\Ji+diq{D, 2)Ji{d, q) 



{/-l)g ] [~T^] pT"'SD-iSd-iJM/giD,2)Ji{d,q) . (A.9) 
Now we use mcoj/k = r^^^/4pg and 

mta^/2 (moDj^''' 1 1 Nt - I 



to write 



and 



KAN,d,q) = ^^iw/.(D,2)7K^,^) ^^^^-^,,,^ ^^j -f^ (A.11) 



is-^A?^ ^ s-^w^ , JM/q{d,q)Ji{d,q) / \ 

^:/(A^, J, ^) - i«:i(A^, d, q)- , , . . , , . — — — . (A. 12) 

Ji+dlq{d,q)J\{d,q) \{N - \)g 



[1] I. E. Mandelstam and I. Tamm, J. Phys. (USSR) 9, 249 (1945). 



32 



[2] C. W. Helstrom, Quantum detection and estimation theory, vol. 123 of Mathematics in science and engineering 

(Academic Press, New York, 1976), 1st ed. 
[3] A. S. Holevo, Probabilistic and statistical aspects of quantum theory, vol. 1 of North-Holland series in statistics 

and Probability theory (North-Holland, Amsterdam, 1982), 1st ed. 
[4] S. L. Braunstein and C. M. Caves, Phys. Rev. Lett 72, 3439 (1994). 
[5] S. L. Braunstein, C. M. Caves, and G. J. Milbum, Ann. Phys. 247, 135 (1996). 
[6] V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. Lett. 96, 010401 (2006). 
[7] S. Boixo, S. T. Flammia, C. M. Caves, and J. Geremia, Phys. Rev. Lett. 98, 090401 (2007). 
[8] S. Boixo, A. Datta, S. T. Hammia, A. Shaji, E. Bagan, and C. M. Caves, Phys. Rev. A 77, 012317 (2008). 
[9] S. Gleyzes, S. Kuhr, C. Guerlin, J. Bemu, S. DelegUse, U. B. HofF, M. Brune, J.-M. Raimond, and S. Haroche, 

Nature 446, 297 (2007), ISSN 0028-0836. 
[10] J. J. BolUnger, W. M. Itano, D. J. Wineland, and D. J. Heinzen, Phys. Rev. A 54, R4649 (1996). 
[1 1] S. F. Huelga, C. Macchiavello, T. PelUzzari, A. K. Ekert, M. B. Plenio, and J. I. Cirac, Phys. Rev. Lett. 79, 3865 

(1997). 

[12] A. Luis, Phys. Lett. A 329, 8 (2004). 

[13] J. Beltran and A. Luis, Phys. Rev. A 72, 045801 (2005). 

[14] A. Luis, Phys. Rev. A 76, 035801 (2007). 

[15] A. M. Rey, L. Jiang, and M. D. Lukin, Phys. Rev. A 76, 053617 (2007). 

[16] S. Choi and B. Sundaram, Phys. Rev. A 77, 053613 (2008). 

[17] M. J. Woolley, G. J. Milbum, and C. M. Caves, eprint, 0804.4540 (2008). 

[18] S. Boixo, A. Datta, M. J. Davis, S. T. Flammia, A. Shaji, and C. M. Caves, Phys. Rev. Lett. 101, 040403 (2008). 
[19] A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 

[20] A. L. Fetter, in Bose-Einstein condensation in atomic gases, edited by M. Inguscio, S. Stringari, and C. Wieman 
(lOS press, Amsterdam, 1999), vol. 140 of International School of Physics Enrico Fermi, p. 201. 

[21] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Reviews of Modem Physics 71, 463 (1999). 

[22] J. Rogel-Salazar, S. Choi, G. H. C. New, and K. Bumett, Joumal of Optics B: Quantum and Semiclassical Optics 
6, R33 (2004), ISSN 1464-4266. 

[23] J. E. Williams, Ph.D. thesis. University of Colorado (1999). 

[24] M. R. Matthews, D. S. Hall, D. S. Jin, J. R. Ensher, C. E. Wieman, E. A. Comell, F. Dalfovo, C. Minniti, and 

S. Stringari, Phys. Rev. Lett. 81, 243 (1998). 
[25] D. S. Hall, M. R. Matthews, J. R. Ensher, C. E. Wieman, and E. A. Comell, Phys. Rev. Lett. 81, 1539 (1998). 
[26] D. S. Hall, M. R. Matthews, C. E. Wieman, and E. A. Comell, Phys. Rev. Lett. 81, 1543 (1998). 
[27] K. M. Mertes, J. W. Merrill, R. Carretero-Gonzalez, D. J. Frantzeskakis, P. G. Kevrekidis, and D. S. Hall, Phys. 

Rev. Lett. 99, 190402 (2007). 
[28] R. P Anderson, C. Ticknor, A. I. Sidorov, and B. V. Hall, Phys. Rev. A 80, 023603 (2009). 
[29] A. Shaji and C. M. Caves, Phys. Rev. A 76, 0321 1 1 (2007). 



